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Abstract 

We analyze the potentiality of the next generation of e^e~ linear colliders 
to search for large extra dimensions via the production of fermions pairs in 
association with Kaluza-Klein gravitons (G), i.e. e~^e~ ffG. This process 
leads to a final state exhibiting a significant amount of missing energy in 
addition to acoplanar lepton or jet pairs. We study in detail this reaction 
using the full tree level contributions due to the graviton emission and the 
standard model backgrounds. After choosing the cuts to enhance the signal, 
we show that a linear collider with a center-of-mass energy of 500 GeV will 
be able to probe quantum gravity scales from 0.96 (0.86) up to 4.1 (3.3) TeV 
at 2 (5)(T level, depending on the number of extra dimensions. 



I. INTRODUCTION 

Recently there has been a great interest in the possibility that the scale of quantum 
gravity is of the order of the electroweak scale instead of the Planck scale Mpi ~ 10^^ 
GeV. A simple argument based on the Gauss's law in arbitrary dimensions shows that the 
Planck scale is related to the radius of compactification (R) of the n extra dimensions by 

1 



Ml ~ i?"M^+2 , (1) 

where Ms is the (4 + ra)— dimensional fundamental Planck scale or the string scale. Thus, 
the largeness of the 4— dimensional Planck scale Mpi (or smallness of the Newton's constant) 
can be attributed to the existence of large extra dimensions of volume i?". If Ms ~ 0{1 
TeV), this scenario solves the original gauge hierarchy problem between the weak scale and 
the fundamental Planck scale, and leads to a rich low energy phenomenology []T],3. The 
n = 1 case corresponds to ~ 10*^ km for Ms = 1 TeV, which is ruled out by observations 
of planetary motions. On the other hand, the gravitational force has been shown to 
hold at separations ranging down to 218 fim which corresponds to the bound Ms ^ 3.5 TeV 
forn = 2 [|. 

Fields propagating in compactified extra dimensions give rise to towers of Kaluza-Klein 
(KK) states separated in mass by 0{1/R) 0]. In order to evade strong constraints on 
theories with large extra dimensions from electroweak precision measurements, the Standard 
Model (SM) fields are assumed to live on a 4-dimensional hyper-surface, and only gravity 
propagates in the extra dimensions. This assumption is based on new developments in string 
theory [§-0] . If gravity becomes strong at the TeV scale, KK gravitons should play a role in 
high-energy particle collisions, either being radiated or as a virtual exchange stateQ P,|9Hl6 



The KK gravitational modes contain spin-2, spin-1, and spin-0 excitations, however, the 
spin-1 modes do not couple to the SM particles while the scalar modes couple to the trace 
of the energy-momentum tensor, therefore not interacting with massless particles. 

In this work we study the potentiality of e~^e~ linear colliders (LC) to probe extra 
dimensions through the clean and easy-to-reconstruct process 

e+e--.//^, (2) 

where the missing energy is due to KK graviton radiation and / can be either a muon or 



In the case where brane torsion is considered, which is not mandatory in this model, one can 
obtain strong bounds for the string scale |^. 
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a quark. Since the fermion masses are negligible, only the spin-2 KK modes are relevant. 
This process contains not only the Z*G associated production but also its interference with 
the 7*(j contribution, generalizing the analyses of Ref. Moreover, we not only apply 
realistic cuts and include detector resolution effects, but also take into account irreducible 
backgrounds. 

In the framework of the SM, a final state containing lepton pairs plus missing momentum 
is due to the production of Uuu via ZZ / 7 or WW intermediate states. In the case of jet pairs 
and missing momentum, there is a large additional contribution due to the production of 
qq'lv where the extra charged lepton is lost in the beam pipe. The SM also gives rise to two 
reducible backgrounds via e^e~ ff or 77 ff where the momenta of the fermion pair 
is mismeasured. However, these two backgrounds turn out to be negligible after applying 
our cuts. 

The real KK graviton emission does not interfere with the SM processes, and conse- 
quently, the signal cross section is proportional to M^"~^, with the proportionality constant 
depending on the number of extra dimensions due the sum over the KK modes. The spin-2 
KK graviton radiation gives raise to a very characteristic spectrum of missing energy, and 
we exploit this feature together with other kinematical distributions to enhance the signal 
over the SM backgrounds. 

This paper is organized as follows. In Section |I|, we describe the techniques used to 
evaluate the relevant cross sections, while Section |T| contains the main characteristics of 
the signal and backgrounds and the cuts chosen to enhance the KK graviton radiation 



process. In Section IV, we present our results and conclusions. 



II. CALCULATIONAL TOOLS 

We are considering the production of fermion pairs accompanied by large energy and 
momentum imbalances which can be generated by the emission of undetectable KK gravi- 
tons. In the case of lepton pairs we analyze only the final state involving muons since the 
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inclusion of taus will present just a small gain in the limits due to the loss in the detection 
efficiency. On the other hand, the process involving electrons is described by an other class 
of Feynman diagrams besides the ones computed here, thus deserving an analysis by itself. 

The signal and backgrounds were simulated with full tree level matrix elements and, for 
the hadronic case, we took into account all quark contributions with the exception of the 
top quark. The spin-2 KK graviton emission is described by 14 Feynman diagrams, see Fig. 
m, where the KK graviton is attached to each of the SM fields and vertices appearing in the 
SM process e~^e~ ff. Our notations and Feynman rules for the KK graviton interactions 
are the ones in Ref. |T^, and we have calculated the signal cross section analytically using 



FORM [|T8|. On the other hand, the background matrix elements for e^e^ f f^^ and qq'ii' 
were generated using the package Madgraph [|19[. For the background process e^e" qq'iu 
special care is required in the phase space integration due to a collinear divergence that 
appears when the final lepton is an electron or a positron. In this case we employed the 
prescription for the phase space given in Ref. ||20|| . 



In our analyses, the kinematical region of the // pair is such that it is possible that they 
originate from the on-shell production of Z's, therefore, we must regulate the Z exchange 
diagrams in such a way that we do not spoil the gauge invariance of the scattering amplitude. 
We chose to regulate the Z propagator by introducing finite width effects. In principle, this 
presents an apparent problem as the matrix elements are not gauge invariant when we simply 
add an imaginary part in the the Z propagator. Though a formal prescription should be 
able to take care of this problem, we resorted to an approximation that is extremely reliable 



in such cases [0,^. The prescription is to multiply an overall factor of the form 

(3) 



to the summed squared matrix element and not to include the finite width at the matrix 
element level. Here, s stands for the total center-of-mass energy while rriff is the invariant 
mass of the final state fermion pair. This factor introduces a Breit-Wigner resonance for 
the Z boson and for regions far away from the Z resonance it is essentially of the order of 



unity. 

In order to check our signal matrix element, we verified that it is gauge invariant. Writing 
the total Z (7) amplitude as A^"'^e*^(A;), where eapik) is the polarization vector of a KK 
graviton of momentum k, we explicitly showed that A^^^/cq independently vanishes for both 
Z and 7 exchange diagrams. Furthermore, we used the gauge invariance to simplify the 
expression for the cross section, though it remained huge. 

We included in our analyses the energy losses due to the emission of photons off the 



initial state, which we treated in the structure function formalism The differential 

cross section is then given by 

da = j dXidX2 fe\e{Xl, y/s) fe\eix2, ^/s) da{s = X1X2S) , (4) 

where a{s) is the cross section in the absence of initial state radiation, and 

/e|e(x, y/s) = (3 



(5) 



with (3 = (log ^ — 1^ , is the leading-log re-summed effective distribution function. 
Note that we did not include beamstrahlung, which is expected to further smear out the 
peak in the e~^e~ luminosity at s = s, since it depends on details of accelerator design. 

For strongly interacting final states, we simulated the resolution of the hadronic calorime- 
ter by smearing the energies (but not directions) of all final state partons with a Gaussian 
error given by 



6E 



had V E 



0.30 

0.01 , (6) 



where E is given in GeV. 



III. SIGNAL AND BACKGROUND PROPERTIES 

We started our analyses imposing the following set of acceptance cuts: (CI) We required 
that the events present a missing transverse momentum^ > 10 GeV; (C2) the muons or jets 
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should have a transverse momentum px > 5 GeV; (C3) we also required that the muons/jets 
are observed in the region | cos6\ < 0.98 , where 6 is the muon or jet polar angle; (C4) the 
jets (muons) are required to be separated by AR > 0.4, where (AR)^ = {Ai])^ + (A0)^, 
with T] being the pseudo-rapidity and (p the azimuthal angle. 

In Table |, we display the total signal cross section after these initial cuts for a center- 
of-mass energy of 500 GeV, assuming Ms = 1 TeV and the number of extra dimensions 
n = 2-7. We can see from this table that the signal cross section drops quickly as we increase 
n, as expected. After the acceptance cuts (C1)-(C4), we found that the cross section of the 
SM background e+e~ f f^i^ is cr'^^fj^ = 73.6 fb in the muonic case and ajj^ = 285.8 fb 
in the hadronic case. For the hadronic final state, the presence of a collinear singularity 
leads to a large cross section for the process e~^e~ qq'ev where the is lost in the beam 
pipe. After cuts (C1)-(C4) and requiring | cos6'e+(e-)| > 0.98 in order to miss the into the 
beam pipe, the cross section for this reaction is a'j'f^ = 3276 fb, making this the dominant 
background in the hadronic case. Consequently, we needed to introduce further cuts to take 
care of this background, rendering the analyses of the hadronic case quite different from the 
muonic one. We finally checked that when a muon or tau is lost into the beam pipe instead 
of an e^, this reaction does not give a significant contribution to the total background after 
our acceptance cuts due to the absence of collinear divergences. 

The dramatic drop of the signal cross section for large values of n compels us to refine 
our analyses by studying kinematical distributions in order to determine further cuts to 
enhance the signal. In Fig. we display the missing invariant mass (Mmiss) distributions of 
the SM backgrounds and KK graviton emission signal for the muonic case after imposing 
the acceptance cuts and taking Ms = 1 TeV and n = 3. Although they are experimentally 
indistinguishable, we show the SM background distributions classified by the final neutrino 
flavor in order to see in detail their behaviors. The electron and tau neutrino flavors present 
a peak near the Z mass, since they originate from the ZZ contribution to this process, while 
the distribution due to the muon neutrino is broader, because of the contribution of two W 
production to this final state. Moreover, the electron neutrino background mimics the signal 
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distribution due to W fusion process. Of course, we added up all backgrounds to obtain the 
final results. Notice that the dips in the distributions at high Mmiss values are due to the 
acceptance cuts. Clearly, Mmiss will be an important variable to reduce the SM background, 
and consequently we required (C5) Mmiss > 320 GeV. 

Due to the KK graviton emission the final state jets and muons are not expected to 
be back-to-back. Fig. |^ contains the distribution of the cosine of the angle between the 
final state muons (cos6'^^) after imposing cuts (C1)-(C5) for Ms = 1 TeV and n = 3. We 
can see that the KK graviton signal slightly prefers the region where the two muons are 
close together while the background receives a large contribution from muons in opposite 
hemispheres. Therefore, we further demanded that (C6) cos6^^ > 0. 

We display in Table || the total cross section for the muonic signal after cuts (C1)-(C6) 
as a function of the number of large extra dimensions for Ms = 1 TeV. After these cuts, the 
SM background is reduced to ct|^^ = 3.24 fb, a reduction by a factor of more than 20 while 
the signal is reduced by a factor of less then 3. Therefore, these cuts enhance considerably 
the signal and extend the attainable bounds for the number of extra dimensions up to six. 

In the hadronic case the main background is e~^e~ — > qq'ev. This background is mainly 
due to Wev production with the W further decaying in qq' . Thus a cut in the invariant 
mass of the jet pair near the W mass can reduce substantially this background. However, 
the signal cross section peaks at the Z mass, as we can see from Fig. ^. Nevertheless, we 
found a cut in the jet pair invariant mass that suppresses the background while exhibiting 
a good efficiency for the signal |2|]: (C7) Mjj < 35 GeV or M^j > 85 GeV. 

In the background qq'eu, a large amount of the missing energy is carried by the 
escaping through the beam pipe. To illustrate this fact, we present in Fig. | the distribution 
of the cosine of the polar angle of the missing momentum after applying the acceptance cuts 
(C1)-(C4), which suggested the introduction of the following cut: (C8) | cos 9miss\ < 0.8. 

After these cuts the qq'uu total cross section is comparable with qq'eu; see Fig. ^ which 
was obtained applying cuts (C1)-(C4) and (C7)-(C8) for ^ = 500 GeV, Ms = 1 TeV, and 
n = 3. The qq'vv reaction presents a peak in the missing invariant mass {Mmiss) distribution 
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near the Z mass. In order to reduce this background we imposed the conservative cut (C9) 
M„,,, > 200 GeV. 

The effects of cuts (C1)-(C4) and (C7)-(C9) in the signal cross section are shown in Table 
^ where we can see that the total background is reduced from o^-^^ = 3562 fb to 284.7 fb, 
while the signal is reduced by a factor of two at most, resulting in significant enhancement 
of signal over background. The effects of cuts on the signal are larger for smaller number of 
dimensions since the cross section for each KK mode goes as M^'^~^mjir^, with being 
the KK mode mass. However, for n > 5 the ratio S/B is lower than 0.1, which means that 
we need a precise background estimation. In order to have a better S/ B ratio we imposed 
more stringent cuts. We found out that requiring (CIO) a harder missing invariant mass 
of Mmiss > 300 GeV; and (Cll) an acoplanarity cut of cosOjj > 0.8, produces the results 
shown in the last line of Table ||. As we can see, even for n = 6 we have S/B > 0.1. 

IV. RESULTS AND DISCUSSIONS 

In this work we studied the potential of e^e~ colliders to probe the quantum gravity scale 
Ms via the KK graviton emission associated with two fermions (muons and quarks). We 
considered a LC with a center-of-mass energy ^/s = 500 GeV and three different integrated 
luminosities C = 50, 200, and 500 fb~^. We derived the constraints on Ms assuming that 
no deviation from the SM predictions was observed. For a given integrated luminosity C, 
the statistical significance of the signal is 

-signal 



a" 



back 



v^, (7) 



where a''"''^^ is the total SM background cross section and (j'^*^""' the total signal cross section; 
see Sec. |T|. Since the quantum gravity signal does not interfere with the SM backgrounds, 
the KK graviton emission cross section is proportional to M^"~^ and we can write the signal 
cross section as 



a^'^-'^\Ms, V~s) = ^CT^^^"'^'(1 TeV, V~s) , (8) 



where Ms is given in TeV. Notice that cr^^^naz depends upon the number of extra dimensions 
n due to the sum over the KK modes. Therefore, one can obtain the 2 (5)(t bounds on Ms 
from d^) and (H) as 



VCa"^f""'(l TeV, yi)' 



l/(n+2) 



Ms < \ ' . (9) 



We present in Table |T| the 2a and 5a attainable bounds on Ms for several choices of 
n, taking into account the muonic, hadronic, and combined channels and using Eq. (P) and 
the results in Table |I|. The combined limits were obtained requiring S/B > 0.1, except for 
n = 7, and we used the formalism, i.e., the sum of muonic and hadronic should result 



in a Ax^ — 4 (25) to be consistent with 2 {5)a limits as prescribed in p4[. Notice that the 



hadronic bounds are slightly better than the muonic ones for n < 5 and that for small values 
of n, the set of cuts A gives a better significance even though S/ B is bigger for the B cuts. 
Our bounds are comparable with the ones derived from the analysis of the associated KK 



graviton production with a Z in , however, our simulations of the signal and backgrounds 
are more detailed, e.g. we take into account the 7*/Z* interference, more realistic cuts, and 
reducible backgrounds. Moreover, we were able to extend the bounds to all values of n in 
contrast with the results in [|r^ which are valid only for n = 2. It is important to notice 
that for n = 7 the ratio S/ B 0.05 which means that the bounds presented here should be 
taken with a pitch of salt since we need a more careful study of the backgrounds in order to 
take seriously the n = 7 limits into account. 



We can also see from Table |ITT| that a 500 GeV LC will be able to probe the quantum 
gravity scale Ms above the new gravitational direct experimental limit Ms > 3.5 TeV for 
n = 2 provided its integrated luminosity is larger than 200 fb~^. Present collider limits 
are, in general, less stringent than this one. For instance, graviton direct production at 
CERN 200 GeV LEP gives 95% C.L. limits of Ms > 1.02 - 1.25 TeV for n = 2 E^, and 



searches for virtual graviton effects lead to 95% C.L. limits of Ms > 0.75 - 1.3 TeV for any 



number of extra dimensions depending on the LEP experiment EH]. The Fermilab Tevatron 



D0 experiment presented a 95% C.L. bound of Ms > 1.37 TeV for n = 2 based on searches 



for virtual graviton effects on dielectron or diphoton systems p6 



We should compare our limits to the ones obtained from alternative signatures in 500 
GeV e^e^ LC. The most significant bounds come from 76* production giving a 5cr limit 



of Ms > 3.66 TeV and gauge boson pair production {VV) giving a 2a limit of Ms > 
2.8 (3.0) (3.2) TeV ioi V = W (Z) (7) []T1|. These results are comparable to the ones 



presented in Table |III|, however we expect that our signature can give direct information 
about the graviton spin through the study of the angular distributions of the final particles. 

In brief, we showed in detail that 500 GeV LC with an integrated luminosity of 500 
fb^^ will be able to exclude Ms up to 4.1 (1.1) TeV for n = 2 (7). Although our results 
are well above the actual experimental limits from LEP and Tevatron, they are a factor 
of 2 less stringent than the expected ones from the CERN Large Hadron Collider (LHC) 



15 , however the LC leads to cleaner and easier to reconstruct events. Moreover, at the 



LHC there is an ambiguity on how to unitarize the cross sections since at very high parton- 
parton center-of-mass energies the subprocesses involving KK gravitons violate unitarity. 
The signal studied here are free from this ambiguity as we have the direct production of 
the graviton. In addition it might be possible to probe the spin of the graviton looking at 
angular distributions PB . 
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TABLES 



n 


2 


3 


4 


5 


G 


7 




55.1 


17.2 


6.08 


2.27 


0.888 


0.357 


(fb) 


723. 


203. 


64.4 


21.9 


7.82 


2.86 



TABLE I. Total signal cross section in fb for the muonic channel ("T^*^" ) and hadronic channel 
(cr^*^"'"') for different number of extra dimensions, using -y/s = 500 GeV and Ms = 1 TeV after 
applying the acceptance cuts (C1)-(C4), as explained in the text. For comparison, the total cross 
sections for the SM backgrounds are = 73.6 fb and a''%'' = 3562 fb. 



n 


2 


3 


4 


5 


6 


7 


(fb) 


18.7 


7.46 


3.05 


1.27 


0.537 


0.230 


(fb) - cut A 


387. 


123. 


40.7 


14.1 


5.06 


1.87 




30.5 


12.4 


5.16 


2.18 


0.931 


0.400 



TABLE II. Total signal cross sections in fb for the muonic and hadronic channels as a function 
of the number of extra dimensions, assuming -s/s = 500 GeV and Ms = 1 TeV. In the muonic 
case, we applied the cuts (C1)-(C6). In the hadronic case, set of cuts A stands for cuts (C1)-(C4) 
and (C7)-(C9) while set of cuts B for (C1)-(C4) and (ClO)-(Cll). The SM cross sections are 
= 3.24 fb and = 285. (8.13) fb for cuts A (B). 
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n 


2 3 4 5 6 7 


muonic 


50 fb-i 
200 fb^i 
500 fb-^ 


2.46(1.96) 1.71(1.42) 1.35(1.16) 1.14(1.00) 1.01(0.90) 0.92(0.83) 
2.93(2.33) 1.97(1.64) 1.51(1.30) 1.26(1.10) 1.10(0.98) 0.99(0.89) 
3 28l'2 61'l 2 15('l 79) 1 63('l 40) I M(\ \H) 1 16fl 04) I 04(0 94") 

KJt^Kjy^tKJj^j ^•xt_/vi.ity/ x.wtj\-L«T:\j/ x.(j^\i.iu/ x.xvj\-L«v^^/ x.w^\V7«tyT:/ 


hadronic 


cut A 


50 fb-i 
200 fb-i 
500 fb-i 


3.00(2.39) 1.91(1.59) 1.43(1.23) 1.17(1.02) 1.01(0.90) 0.90(0.81) 
3.57(2.84) 2.20(1.83) 1.60(1.38) 1.29(1.13) 1.10(0.98) 0.97(0.88) 
4 oof's! 8^1 '?4^(2 0^) 1 73('1 49'> 1 38n 21 1 ^(^(^ 04) 1 01(0 Ql) 

t:.vJUI(J.XOJ ^.T:i\^.VJX/ X.it_l\i.^t7/ X.(JO\-L.^X/ X.XVJVi.vJ^/ X.U^yv^.ty^/ 


cut B 


50 fb-i 
900 fh~l 
500 fb-i 


2.48(1.97) 1.73(1.44) 1.36(1.17) 1.15(1.01) 1.02(0.91) 0.93(0.84) 
2.95(2.35) 1.99(1.65) 1.53(1.31) 1.27(1.12) 1.11(0.99) 1.00(0.90) 
3.31(2.63) 2.18(1.81) 1.65(1.42) 1.36(1.19) 1.18(1.05) 1.05(0.95) 


combined 


50 fb-i 
200 fb^i 
500 fb"^ 


3.07 (2.44) 1.97 (1.64) 1.48 (1.27) 1.20 (1.06) 1.06 (0.94) 0.96 (0.86) 
3.65 (2.91) 2.26 (1.88) 1.66 (1.42) 1.33 (1.17) 1.15 (1.03) 1.03 (0.93) 
4.10 (3.26) 2.48 (2.06) 1.79 (1.54) 1.42 (1.25) 1.22 (1.09) 1.09 (0.98) 



TABLE III. 2 (5)(T limits in TeV on the quantum gravity scale Ms in TeV as a function of 
the number of extra dimensions n for a 500 GeV LC with luminosities C = 50, 200, and 500 fb^^. 
We present the results for the muonic, hadronic with the two cut selections as in Table ||, and 
combined channels. For the latter we required S/B > 0.1, except for n = 7. 
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FIG. 2. Missing invariant mass (Mmiss) spectrum originated from the KK graviton radiation 
(solid line) and the SM contributions to the muonic channel divided in the neutrino flavors: 
(dot-dashed), (dotted), and Ur (dashed). We assumed ^/s = 500 GeV, Ms = 1 TeV, n = 3, and 
applied the acceptance cuts (C1)-(C4) described in the text. 




FIG. 3. Distribution of the cosine of the angle between the final state muons originated from 
the KK graviton radiation and the SM contributions as in Fig. ^ for -y/s = 500 GeV. We assumed 
Ms = 1 TeV, n = 3, and applied the cuts (C1)-(C5) described in the text. The Vr distribution is 
not displayed since it is too small. 
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FIG. 4. Dijet invariant mass (Mjj) spectrum originating from the KK graviton radiation (solid 
line) and the SM contributions to the hadronic case qq'eu (dashed) and qqi'i' (dotted) for -^5 = 500 
GeV. We assumed Ms = 1 TeV, n = 3, and applied the acceptance cuts (C1)-(C4). 




FIG. 5. Missing momentum polar angle {cos Omiss) spectrum coming from KK graviton radia- 
tion and the SM contributions to the hadronic case. We assumed ^/s = 500 GeV, Ms = 1 TeV, 
n = 3, and applied the acceptance cuts (C1)"(C4). The lines follow the conventions as in Fig. ^. 
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FIG. 6. Missing invariant mass {Mmiss) spectrum originating from KK graviton radiation and 
the SM contributions to the hadronic case. We assumed ^/s = 500 GeV, Ms = 1 TeV, n = 3, and 
apphed the acceptance cuts (C1)-(C4) plus the kinematical cuts (C7)-(C8). The hues fohow the 
convention as in Fig. ^ 
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